Melting of the stripe phases in the t-t'-U Hubbard model 
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Abstract 

We investigate melting of stripe phases in the overdoped regime x ^ 0.3 of the two-dimensional t-t'-U Hubbard 
model, using a spin rotation invariant form of the slave boson representation. We show that the spin and charge order 
disappear simultaneously, and discuss a mechanism stabilizing bond-centered and site-centered stripe structures. 
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It is now well established that the doped cuprates 
show many highly unusual properties both in normal 
and superconducting state. Among them, stripe phase, 
discovered in theory [1] and confirmed by experiment 
[2], attracted a lot of interest. Instead of moving in- 
dependently, the holes introduced to an antiferromag- 
netic (AF) Mott insulator self-organize either on site- 
centered (SC) nonmagnetic domain walls (DW) sep- 
arating AF spin domains, or on bond-centered (BC) 
DW made out of pairs of ferromagnetic spins [2] . Such 
a tendency towards phase separation is fascinating, 
and offers a framework for interpreting a broad class 
of experiments, including the pseudogap at the Fermi 
energy observed in the angle-resolved photoemission 
(ARPES) spectra of La 2 _ I Sr ;c Cu04 (LSCO) for the 
entire underdoped regime (0.05 ^ x ^ 0.125), repro- 
duced within the t-t'-U Hubbard model [3]. Therefore, 
we argue that this model is sufficient to investigate 
generic features of stripe phases. 

Two main scenarios for a driving mechanism of the 
stripe phase have been proposed [4]. In the first one 
stripes arise from a Fermi surface instability with the 
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spin driven transition [1]; then spin and charge or- 
der simultaneously, or charge order follows spin order. 
The second scenario comes from Coulomb-frustrated 
phase separation suggesting that stripe formation is 
commonly charge driven, and the charge order sets in 
first when the temperature is lowered. However, slave 
boson studies of the two-dimensional (2D) t-t'-U Hub- 
bard model showed that the spin susceptibility diverges 
while the charge susceptibility does not [5] , so the mi- 
croscopic origin of the stripe instability is unclear. 

We investigate the mechanism leading to phase sep- 
aration and the melting of vertical BC and SC stripe 
phases in the overdoped regime (x ^ 0.3, where x = 
1 — n and n is an average electron density per site) 
of the 2D t-t'-U Hubbard model. We employ the spin 
rotation invariant slave boson (SB) representation of 
the Hubbard model [6], and perform the calculations 
on larger (up to 144x144) clusters than those studied 
recently [7] . This allows one to obtain unbiased results 
at low temperature T — O.Olt. 

For the model parameters for LSCO: U/t = 12 and 
t' It — — 0.15, we obtain that the most stable SC stripes 
are separated by d — 4(3) lattice spacings at dopings 
0.124 < x sC 0.2 (0.2 sC x < 0.34), respectively. As 
shown in Fig. 1(a), increasing doping stabilizes the SC 
stripes with a single atom in the AF domains. Also for 
the BC stripes the size of the AF domains decreases 
with increasing doping, varying from d = 5 (0.10 ^ 
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Fig. 1. Melting of vertical BC and SC with increasing doping 
x at T — O.Oli: (a) the free energy SF (black line) and inter- 
action (grey line) energy gain in the stripe phases; (b) local 
charge densities Srii relative to their average values; (c) local 
magnetization m^; (d) double occupancies SDi relative to the 
values in the paramagnetic phase (scaled by a factor for the 
d — 2 stripe). In panels (b)-(d) the black (grey) curves corre- 
spond to the strongly (weakly) polarized sites, respectively. 

x 0.13) through d = 4 (0.13 < x < 0.19) and down 
to d = 3 at higher doping, as there is no BC configura- 
tion with d — 2. For both types of stripes, the distance 
between them is locked to four in a sizeable doping 
range above x ~ i, in agreement with neutron scat- 
tering experiment [8] and with theory [3] for LSCO. 

In Fig. 1(a) we show the energy gain of the stripe 
phases with respect to the paramagnetic phase SF. Re- 
markably, the difference in energy between the best SC 
and BC stripes is smaller than both the accuracy of the 
calculations, and the resolution of Fig. 1(a), suggest- 
ing that quantum fluctuations might be important. We 
characterize the melting of stripes by their SB local 
averages: density m = zCr^f), magnetization rrii = 
| (Si) |, and double occupancies Di = {n^nn). 

In the d = 2 SC stripe, reported here for the 
first time, the two 5rii(x) curves are symmetrical in 
Fig. 1(b). In contrast, in the d = 3 BC stripe there 
are two sites with weak magnetic moments per one 
strongly polarized site. We note that, unlike in the 
SC phase, the variation in density is largest on the 
strongly polarized sites in the BC phase. The mag- 
netic moments rm vanish for both types d — 3 stripes 
at the same doping x = 0.375 [Fig. 1(c)], suggesting 
that they originate from the same instability. 

The microscopic mechanism stabilizing the d = 2 SC 
stripes appears to differ markedly from the one stabi- 
lizing the d = 3 ones [9]. For d = 2 [Fig. 1(d)], the re- 
duction of double occupancy is strongest on the mag- 



netic sites, and the corresponding reduction of interac- 
tion energy is larger than the gain of free energy [see 
Fig. 1(a)]. Thus the mechanism leading to the forma- 
tion of the d — 2 stripe is primarily local, making use 
of two complementary effects helping to reduce double 
occupancy: finite magnetization at magnetic sites and 
reduced electron density at nonmagnetic ones. Even 
though such a state looses kinetic energy, the gain in 
the interaction energy overcompensates this loss, sta- 
bilizing this order in a wide doping range x ^ 0.485. 

In contrast, for d — 3 stripes, both contributions 
to the free energy are substantially decreased while 
stripe order starts melting already at x < 0.3 mainly 
by faster removing double occupancies from the stripe 
DW than from AF domains leading to gradually dis- 
apearing magnetic moment upon doping. Therefore, 
both potential and kinetic energy (including the su- 
perexchange) cooperate to stabilize stripes with d > 2. 
In fact, for both d = 3 stripes, the mechanism is dop- 
ing dependent. In the small magnetization regime, the 
interaction energy plays the leading role. However, un- 
der a further decrease of hole density, this gain nearly 
saturates (at x ~ 0.33), and the gain in the kinetic en- 
ergy starts to dominate. Moreover, it is only slightly 
larger for the SC stripe compared to the BC one, and 
therefore it is easily compensated, mainly by the pres- 
ence of finite magnetic moments at BC domain walls. 
As a common feature, the spin and charge order disap- 
pear at the same critical doping. Therefore, in the ab- 
sence of longer ranged Coulomb interaction the charge 
order is always accompanied by the spin order. 

Summarizing, we have investigated the microscopic 
mechanisms responsible for the formation of the ver- 
tical BC and SC stripes in the extended 2D Hubbard 
model. Interestingly, we found that BC and SC stripes 
remain nearly degenerate, and both spin and charge 
order vanish simultaneously when they melt, demon- 
strating a cooperative character of the stripe order. 
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